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SUMMARY 

The  object  of  this  thesis  ia  to  determine  the  effect  of  small 
openings  on  the  atresa  distribution  in  webs  of  beams  subjected  to  pure 
bending  and  to  devl8e  a  rapid  means  for  predicting  the  increase  in  stress 
over  that  of  the  condition  of  no  opening. 

Experimental  evidence  for  a  family  of  specimens  which  covered 
the  usual  structural  practice  he  a  been  plotted  in  the  RESULTS.  A  theor- 
etical solution*  developed  in  the  APPENDIX,  is  included  and  compared  with 
the  experimental  evidence.  From  these  plots  the  authors  conclude  thatt 
the  theoretical  solution  compares  quit©  favorably  with  actual  experiment! 
the  stress  concentration  factor  for  the  edge  of  the  hole  varies  directly 
with  -the  ratio  of  depth  of  opening  to  depth  of  web  and  inversely  with 
the  ratio  of  length  to  depth  of  opening;  and  that  the  stress  concentration 
factor  for  the  edge  of  the  web  varies  directly  with  the  ratio  of  depth 
of  opening  to  depth  of  web,  but  the  variation  with  the  ratio  of  length 
to  depth  of  opening  is  more  complex  showing  a  pronounced  maximum  in  the 
vicinity  of  1.5.  The  fact  that  the  ratio  of  depth  of  opening  to  depth 
of  web  ia  of  little  importance  until  it  reaches  the  value  of  about  1/3 
is  considered  excellent  verification  of  presently  accepted  practice. 

On  the  basis  of  the  foregoing,  it  is  considered  highly  desirable 
that  further  investigation  for  specimens  having  a  ratio  of  length  to  depth 
of  opening  in  the  neighborhood  of  1.5  be  rn.-i.de  with  a  view  to  confirmation. 
It  is  further  considered  desirable  that  a  systematic  survey  of  a  similar 
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family  be  made  to  determine  the  effectiveness  of  reinforcement  at  open- 
ings and  simultaneously  cases  of  complex  bendingj  i.e.,  bending  with  shear. 


-2- 


IHTRODDCTIOa 

In  almost  all  structures  openings  of  some  sort  are  required, 
either  for  access  of  for  passing  piping,  ventilation,  or  electrical  systems 

to  the  various  compartments.  It,  therefore,  behooves  the  structural 
designer  to  make  proper  provision  for  these  requirements  without  adversely 
affecting  the  structure. 


Haturally,  the  circle,  because  of  its  simplicity,  was  the  first 

* 
(1) 
type  of  opening  chosen.  Kir scir  'analysed  the  effect  of  a  circular 

(2) 
opening  in  a  uniform  tension  fieldj  Tuaiv  '   performed  the  analysis  for 

the  circular  opening  in  pure  bending;  Rowland  and  Stevenson w/  extended 

both  the  Kirsch  and  Tusi  analyses  for  finite  width  of  plate  and  introduced 

U)         (5) 
the  analysis  for  complex  bending.  Wolf   and  Heuber   performed  analyses 

(6) 

for  the  elliptical  opening  in  pure  bendingj  Inglis   made  a  similar  analysis 

for  the  elliptical  hole  in  a  uniform  tension  field;  Wolf  also  extended 
his  analysis  to  cover  the  ellipse  in  complex  bending.  Quite  recently 

Greenspan*  ' introduced  an  analysis  for  the  flat  ovaloid  in  a  uniform 

(8) 
tension  field.  In  a  forthcoming  publication  ,  Joseph  and  Brock  develop 

a  theoretical  solution  for  the  flat  ovaloid  in  pure  bending.  This  theBis 

is  devoted  to  the  latter  problem.  The  theoretical  solution  in  the  AP.BK&IX 

follows  the  method  proposed  in  (8).  The  effect  of  finite  depth  is  evolved 

by  an  analysis  of  experimental  data  rather  than  by  a  rigorous  mathematical 

method. 


lumbers  in  brackets  refer  to  the  Bibliography. 


-3- 


It  is  the  purpose  of  thi3  thesis  to  correlate  data  frost  photo- 
elastic  experiments  wito  a  systematic  series  of  models  with  varying  ratios 
of  length  to  depth  of  opening  and  depth  of  opening  to  depth  of  web,  eub- 
jected  to  pure  bending,  in  which  the  openings  are  flat  ovaloidsj  i.e., 
openings  ccaposed  of  rectangles  vrith  semicircles  mounted  at  tne  long- 
itudinal ends.  This  correlation  includes  a  check  on  the  mathematical 
solutions  developed  in  (2),  (3)#  (8)»  and  in  the  APPENDIX. 
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PROCEDURE 

The  models  were  prepared  using  a  nominal  thickness  of  3/3 
inches  clear  Catalin.  All  models  were  cut  from  the  same  sheet  of  ntaterial, 
and  a  standard  tensile  specimen  was  included.  The  tensile  specimen  was 
used  only  to  determine  the  fringe  constant  of  the  particular  sheet  to 
ensure  satisfactory  sensitivity  as  well  as  to  provide  means  for  obtain- 
ing checks  through  the  equations  of  static  equilibrium,  if  necessary. 
These  models  were  machined  using  high-speed  steel  end  mills  for  all  cutting* 

In  an  effort  to  minimise  the  time-edge  effect,  polymerisation, 
a  speoimen  was  tested  immediately  after  machining,  allowing  only  suffic- 
ient time  for  smchining  stresses  tc  dissipate.  A  white  light  polariscope 
was  used  for  this  purpose.  The  specimen  was  fitted  with  end  brackets 
and  suspended  in  an  A  -  frame,  sa  shown  in  Fig.  I,  whioh  subjected  the 
entire  specimen  to  pure  bending.  Suspension  was  accomplished  by  vertical 
links  fitted  at  each  end  with  circular  pins  of  drill  stock  to  assure  as 
nearly  point  support  and  loading  as  was  feasible.  The  A  -  frame  was  loca- 
ted in  the  center  of  the  circular  polariscope  using  a  mercury  vapor  lamp 
with  a  Wratten  #62  Filter.  The  resulting  light  whs  very  nearly  mono- 
chromatic with  a  wave  length  cf  54ol  A°.  The  resulting  fringes  were 
alternately  light  and  dark  giving  excellent  contrast. 

The  specimen  was  first  suspended  with  brackets  and  no  load* 
The  analyser  was  rotated  to  give  a  dark  field.  A  picture  was  taken  in 
this  condition,  hereinafter  referred  to  as  "no  load".  The  negative 
was  then  developed  to  assure  proper  focus  before  load  was  applied. 
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Load  was  then  applied,  and  a  new  picture  token.  The  load  was  previously 
determined  which  would  give  a  maximum  stress  of  approximately  100  pel., 
well  within  the  elastic  range  yet  sufficiently  high  to  give  about  seven  (7) 
fringes.  The  analyser  remained  in  the  same  position  as  in  the  no  load 
condition.  The  negative  was  then  developed  to  assure  satisfactory  focus 
and  contrast  before  the  specimen  was  removed  for  further  machining.  By 
following  this  method,  it  was  possible  to  complete  a  run,  five  (5)  specimens 
-  all  with  the  same  ratio  of  length  to  depth  of  opening,  but  with  varying 
ratios  of  depth  of  opening  to  depth  of  web,  in  about  four  (4)  hours  which 
minimised  the  time-edge  effect. 
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RESULTS 

The  results  of  *he  experimental  data  are  incorporated  in  Figures 
II,  III,  IV,  and  V.  Figure  III  is  a  plot  of  stress  concentration  factor 
at  edge  of  hole  versus  the  ratio  of  length  to  depth  of  opening  for  contours 
of  depth  of  opening  to  depth  of  web.  Figure  II  is  a  plot  of  stress  con- 
centration factor  at  the  edge  of  the  hole  versus  the  ratio  of  depth  of 
hole  to  depth  of  web  for  contours  of  constant  ratios  of  length  to  depth 
of  opening.  On  this  plot  the  various  experimental  data  is  sliown.  Simil- 
arly Figure  IV  and  T  show  respectively  the  same  plots  with  stress  concen- 
tration factor  at  the  edge  of  web  as  the  substituted  variable.  Table  I 
is  a  compilation  of  experimental  data  from  which  the  aforementioned  curves 
have  been  plotted. 
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flunaary  of   experimental  data  obtained  photoelaatically 

Ratio  of  length     Ratio  of  depth     Stress  Concentration  Factor 
to  depth  of        of  opening  to 
opening  depth  of  beam 
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DISCUSSION  OF  RESULTS 

The  curves  shown  in  Figures  II  and  III,  and  IV  and  V  respectively 
represent  the  cross-fairing  of  the  data  incorporated  in  Table  I*  The  results 
of  Figures  II  and  III  are  consistent  with  the  well-verified  facts  that 
a  small  crack  in  the  longitudinal  direction  of  a  benra  causes  a  stress  con- 
centration of  unity  and  a  small  crack  in  the  transverse  direction  of  the 
bean  causes  an  infinite  stress  concentration.  This  character  is  shown 
again  in  the  slope  of  the  curves  of  Figures  III.  It  is  not  surprising, 
either,  that  the  circular  opening  causes  a  greater  stress  concentration 
than  any  of  the  other  flat  ovaloids  (no  ovaloids  transverse  to  the  axis 
of  the  beam  were  tested) •  This  result  is  similar  in  content  to  a  compar- 
ison of  (1)  with  (2)  and  of  (6)  with  U) . 

The  curves  of  Figures  IV  and  7  are  not  so  easily  explained. 
At  the  outset  there  appears  to  be  no  rational  reason  for  the  apparent 
maximum  stress  concentration  factor  at  the  edge  of  the  beam  for  a  ratio 
of  length  to  depth  of  opening  of  1.5.  The  almost  perfect  agreeaent  of 
the  experimental  results  with  the  theoretical  solution  for  the  value  of 
the  ratio  of  depth  of  opening  to  depth  of  beam  of  sero  ana  for  all  values 
of  suoh  ratio  for  the  circular  opening  lend  confidence  in  the  remaining 
experimental  work.  On  this  basis  the  authors  feel  safe  in  predicting  that 
the  maximum  augmentation  of  stress  caused  by  a  small  opening  of  this  type 
will  occur  when  the  ratio  of  length  to  depth  of  opening  is  in  the  vicinity 
of  1.5.  This  same  confidence  enables  the  authors  to  conclude  that  the 
ratio  of  depth  of  opening  to  depth  of  beam  should  lie  between  1/4.  and  1/3 
if  no  material  increase  in  the  stress  concentration  at  the  edge  of  the 
hole  is  desired. 
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CONCLUSIONS 

1.  The  stress  concentration  factor  at  the  edge  of  the  hole  varies 
inversely  with  the  ratio  of  length  to  depth  of  opening  and  directly  with 
the  ratio  of  depth  of  opening  to  depth  of  beam, 

2.  The  stress  concentration  factor  at  the  edge  of  the  beam  varies 
directly  with  the  ratio  of  depth  of  opening  to  depth  of  beam,  rut  in  seme 
complex  manner  with  the  ratio  of  length  to  depth  of  hole, 

3.  The  critical  value  of  the  ratio  of  length  to  depth  of  opening 
falls  in  the  vicinity  of  1.5-  The  critical  value  of  the  ratio  of  depth 
of  opening  to  depth  of  beam  lies  in  the  neighborhood  of  1/4  to  1/3 •  The 
critical  value  is  that  value  at  vhich  a  small  increase  in  the  aforementioned 
ratio  -"  s  accompanied  by  a  relatively  large  increase  in  the  stress  concen- 
tration factor  at  the  edge  of  the  hole, 

4.  The  theoretical  solution  developed  in  (8)  and  in  txhe  APVENDIX 
is  satisfactory  for  purposes  of  engineering  design. 

5.  The  photoelastic  method  is  well-suited  for  inexpensive  and  rapid 
determination  of  stress  concentration  factors. 
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RECOMMEMDATXUNS 

1.  Further  experimental  investigations  for  values  of  the  ratio 
of  length  to  depth  of  opening  in  the  vicinity  of  1.5  seem  highly  desir- 
able to  fix  this  critical  value. 

2.  Further  experimental  and  theoretical  investigations  in  this 
nature  to  determine  the  effect  of  reinforcement  of  such  openings  and  the 
nature  of  complex  bending;  that  is,  bending  with  shear,  seem  desirable* 

3.  In  the  absence  of  a  olearly  defined  'ra  load*  picture  in  a  photo- 
elastic  investigation,  it  is  recommended  that  small  but  appreciable  stress 
be  introduced  to  clarify  the  basic  pattern. 
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APPESDIX 
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APPENDIX  A 

DETAILS  OF  PROCEDURE 

Description  of  Models 

Pour  (4)  basic  models  were  prepared.  Each  was  fitted  with  an 
interior  hole  1/4.  inches  in  the  transverse  direction,  but  the  longitud- 
inal extent  were  respectively,  1/4  inches,  3/8  inches,  1/2  inches,  and 
5/8  inches.  The  distance  between  suspension  points  was  maintained  at  11  1/8 
inches;  the  distance  from  suspension  points  to  the  loading  points  was 
maintained  at  5  inches.  Bach  of  these  models  was  successively  diminished 
in  dapth  in  order  to  obtain  the  increasing  values  of  the  ratio  of  depth 
of  opening  to  depth  of  beam.  Only  in  the  cajje  of  the  circular  opening 
was  it  necessary  to  recut  the  interior  hole.  The  particulars  of  the  various 
models  and  the  loads  to  which  they  were  subjected  are  summarized  in  Table  II. 
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TABLE  II 

Sumaary  of  Model  Particulars 

Hole  Dimensions 

Depth  of  Beaa 

b/c 

Load  in  lb*. 

a/b 

1/4"  *  iA« 

2.240 

0.1115 

85.35 

1.0 

1/2    -x  1/2 

2.094 

0.239 

85.35 

1.0 

1/2     x  1/2 

1.372 

0.^64 

46.00 

1.0 

1/2     x  1/2 

1.030 

0.485 

22.00 

1.0 

1/2     x  1/2 

0.833 

0.600 

4.00 

1.0 

1/4    *  3/8 

2.500 

0.100 

85.35 

1.5 

1/4    x  3/3 

1.2500 

0.200 

20.00 

1.5 

1/4     x  3/8 

0.833 

0.300 

9.00 

1.5 

1/4    ^  3/3 

0.625 

0.400 

4.0 

1.5 

1/4    x  3/8 

0.500 

0.500 

1.50 

1.5 

1/4    x  1/2 

2.500 

0.100 

85.35 

2.0 

1A    x  1/2 

1.250 

0.200 

20.00 

2.0 

1/4     x  1/2 

0.833 

0.300 

9.00 

2.0 

lA    x  1/2 

0.625 

0.400 

4.00 

2.0 

lA     x  1/2 

0.500 

0.500 

1.50 

2.0 

1/4    x  5/8 

2.500 

0.100 

56.90 

2.5 

1/4    x  5/8 

1.250 

0.200 

18.50 

2.5 

1/4  xx  5/8 

0.833 

0.300 

9.50 

2.5 

1/4     x  5/8 

0.625 

0.400 

6.00 

2.5 

J/4    x,  5/8 

0.5C0 

0.500 

■  5.Q0 

2.f5 

\ 
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Methods  of  Analysis 

Standard  photoelastic  methods  were  used,  supplemented  where 
required  with  direct  analytic  methods.  A  transverse  section  was  chosen 
a  minimum  of  two  lengths  removed  from  the  longitudinal  end  of  the  hole. 
The  condition  of  stress  at  this  section  should,  according  to  St.  Venant 
be  removed  from  the  influence  of  the  hole.  This  section  was  designated, 
accordingly,  "pure  bending*,  see  Pigs.  VI,  VII,  VIII,  and  H,  and  was  the 
reference  condition  for  determining  the  stress  concentration  factors. 

Examination  of  the  photographs,  Figs.  X,  XI,  XII,  and  XIII, 
revealed  that  the  maximum  stress  at  the  edge  of  the  hole  occured  at  its 
vertical  axis  of  symmetry,  designated  "diameter";  whereas  the  maximum 
stress  at  the  edge  of  the  beam  occured  vertically  above  the  longitudinal 
end  of  the  opening,  designated  "end  of  hole."  It  was  therefore  necessary 
to  determine  the  fringe  order  existing  at  the  free  boundaries  in  these 
two  locations.  Since  a  fringe  is  the  locus  of  points  of  constant  differ- 
ence between  principal  stresses,  one  of  which  must  be  sero  at  a  free  boundary, 
no  other  experimental  information  was  necessary  for  the  evaluation  of  actual 
stress  existing  at  these  free  boundaries.  Fortunately,  the  "exact,"  math- 
ematical, stress  analysis  for  pure  bending  also  reveals  that  one  principal 
stress  is  also  zero  everywhere,  which  fact  obviated  any  additional  inform- 
ation for  the  pure  bending  condition  as  well.  Thus  the  desired  stress 
concentration  factor  was  merely  the  ratio  of  the  fringe  order  actually 
measured  at  a  given  point  to  the  fringe  order  determined  at  the  corres- 
ponding depth  in  the  pure  bending  condition.  Both  fringe  orders  referred 


i 


Hef.  (9).  p.  3! 
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to  above  must  be  the  actual  value  deterrained  by  the  photograph  corrected 
for  any  residual  stress  in  the  *no  load"  condition. 

This  correction  is  shown  in  Figures  VI,  VII,  VIII,  and  IX  as 
the  difference  between  the  dotted  and  solid  lines.  The  case  of  pore 
bending  was  obtained  by  plotting  integral  fringe  orders,  solid  lines  in 
the  photographs*  at  the  proper  fraction  of  depth  from  the  bottom.  The 
upper  and  lower  halves  of  the  bean  were  handled  separately,  since,  in  pure 
bending,  one  half  would  be  in  tension  and  the  other  in  compression,  with 
zero  stress  at  the  middle  surface.  The  best  straight  line,  guided  here 
by  the  theoretical  linear  distribution  of  bending  stresses,  was  passed 
through  the  experimental  points  in  each  half  of  the  beam.  It  will  be 
observed  in  Figures  VI,  VII,  VIII,  and  IX  that  in  no  care  did  these  lines 
pass  through  zero  at  the  middle  surface.  Shifts  of  both  halves  of  these 
lines  affected  to  make  the  lines  pass  through  this  zero  point  and  to  inter- 
cept the  same  value  at  the  edges  of  the  beam.  It  will  be  noted  from  the 
figures  that  in  most  cases  this  was  almost  a  parallel  shift  of  lines  with 
approximately  the  same  values  for  both  halves  of  the  beam.  This  amount 
of  shift,  in  all  but  the  shallowest  models,  was  of  the  3a»e  character  as 
the  load  stress  and  less  than  fringe  order  1.  This  represents  the  residual 
stress  in  the  specimen,  which  from  the  density  of  blackness  in  the  "no 
load*  picture,  must  be  very  nearly  constant  over  the  section.  It  may 
be  noted,  however,  that  in  the  shallowest  beamo  there  is  appreciable  resi- 
dual stress,  more  than  a  single  fringe  order  and  of  the  same  character 
to  the  load  stress.  This  condition  is  caused  by  the  bending  effect  of 
the  end   brackets  and  admits  a  more  straightforward  correction  procedure. 
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The  "diameter"  and  "end  of  hole"  curves  have  this  same  correc- 
tion factor  applied  to  them  in  all  cases  except  for  the  circular  opening, 
Figure  VI  and  X.  This  procedure  was  followed  since  the  density  of 
blackness  on  the  "no  load"  picture  appeared  to  be  the  seme  throughout 
these  specimens  except  for  insignificant  surface  blemishes.  The  circular 
openings  retained  machining  stresses  for  a  considerable  period  of  time. 
The  fringes  representing  these  machining  stresses  were  concentric  rings 
which  permitted  a  plot  of  fringe  order  versus  fraction  of  depth  from  the 
bottom,  3ee  Figure  ¥1.  These  values  were  then  added  algebraically  to  the 
experimentally  determined  values,  previously  plotted.  The  character 
of  the  machining  stress  was  easily  determined  by  noting  the  "bulging" 
or  "depressing"  effect  on  the  fringes  in  the  loaded  condition.  For  example 
when  fringes  observed  in  the  compression  half  of  the  beam  were  "depressed" 
toward  the  hole  and  bent  away  from  the  diameter  toward  the  end  of  the  hole, 
thus  increasing  the  fringe  order  at  the  diameter  of  the  hole,  the  machin- 
ing stress  was  obviously  compression  and  should  be  subtracted  from  the 
experimentally  determined  value  to  obtain  the  actual  fringe  order  caused 
by  the  bending  alone.  Conversely,  if,  in  the  compression  half  of  the 
beam,  fringes  wore  "bulged"  away  from  the  edge  of  the  hole,  causing  dips 
on  either  side  of  the  dianeter,  the  machining  stress  was  tension  and 
should  be  added  to  the  experimentally  determined  value  to  obtain  the  actual 
fringe  order  caused  by  the  bending  alone. 

A  polar  plot  to  obtain  the  experimental  fringe  order  at  the 
diameter  of  the  circular  hole  was  used  in  an  effort  to  ascertain  this 
peak  value  with  greater  accuracy.  This  is  shown  in  Figure  ¥1.  The 
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photograph  -ana   tacked  to  a  drawing  board  and  the  center  of  the  hole  deter- 
mined. Radial  lines  were  then  passed  from  the  center  of  the  circle  through 
the  intercepts  of  the  fringes  at  the  circumference.  The  angles  were  then 
transferred  to  the  plot  and  the  correct  value  of  fringe  order  was  plotted 
radially  on  this  line.  Finally,  a  smooth  curve  was  passed  through  the 
points  and  the  value  of  the  fringe  order  determined  on  the  diameter  ex- 
tended. This  gave  the  desired  end  point  of  the  curve. 
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APPE8DU  B 
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APPENDIX  B 

Sample  Calculations 

The  data  are  summarised  in  entirety  in  Table  I  of  RESULTS. 
Figure*  ¥1 ,  ¥11,  ¥111,  and  IX,  referred  to  in  AFPEMBDC  A. .  are  shown  to 
illustrate  the  methods  of  analysis  and  typical  calculations  of  stress  con- 
centration factors.  One  sample  froa  each  series  of  models  Is  included; 
that  is,  one  for  each  ratio  of  length  to  depth  of  opening  and  each  ratio 
of  depth  of  opening  to  depth  of  beam. 
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APPEHDIX  C 

Original  Data 

The  original  data  is  presented  in  the  form  of  photographs  of 
the  "no  load*  and  load  conditions.   Each  page  of  photographs  includes  all 
the  pictures  taken  for  a  given  series  of  aodelst  Figure  X  for  the  ratio 
of  length  to  depth  of  opening  of  1,  Figure  XI  for  the  ratio  of  1.5»  Figure 
XII  for  2,  and  Figure  XIII  for  2.5. 
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b/C  -  0.239     No  Load 


a/b  -  1 

FI6.X 


b/C  -  O.36"     Load 

^^^= 

b/C  -  0.2J9     Load 


b/C  -  0.1115     Lo« 
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FIG.XC 


b/fr  -  0.5     No  Load 


b/C  -  0.5     Load 


b/C  -  0.4     No  Load 


b/C  -  0.4     Load 


b/C  -  0.3     No  Load 


b/C  -  0.3     Load 


b/C  "0.2     No  Load 


1.  f  ii 


fiPMMiJ! 


b/C  -0.2     Load 


b/C  -  0.1       No  Load 


aA  »  1.5 


b/C  «0.1     Load 


-32- 


FlG."3ir 


b/C  -  0.5       No  Load 


b/C  »0.5     Load 


b/C  -0.4       No  Load 


b/G  ■  0.4       Load 


b/C  ■  O.J     No  Load 


b/C  ■  o.)     Load 


b/C  -0.2       No  Load 


b/C  -  0.2     Load 


i 'i&jm  '^  ;  1 


b/C  -  0.1       No  Load 


a/b  -  2 


b/C  =  0.1     Load 


b/C  -  0.5       No  Load 


b/C  ■  0.4     No  Load 


b/C  -0.5       Load 


b/C  -  0.4       Load 


a/b  =  2.5 


FIG.2HI 


b/C  -  O.J       Load 


b/C  «  0.2       Load 


b/C  -0.1       Load 
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FIG."33T«X)  Z- PLANE 


IT  CIRCLE 


FIG. 32(b)         ?  plane 
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APPENDIX   D 

Nomenclature 

(x,y)  =   Cartesian  co-ordinates 

(a  ,P  )  =    Orthogonal  curvilinear  coordinates 

z  =   x  /  iy 

w  =  a   /  ip 

£  =    e* 

h  =    stretch  ratio 

X  =    e1? 

lL  =    angle   in  z-plane  betveen  tangent   to  curve 

p    =    constant  and  x-axis 

s,    t,    R=    real   parameters 

T  «    t  /  Hi  -     G 
'    s 

J,vl  =  potential  functions  of  complex  variable  z 

<p,\|j  =  potential  function  of  complex  variable 

F  =  Airy  stress  function 

a  ,  o  »   a  ,  ag  ■  normal  stresses;  o  is  normal  to  surface  for  *fhich 

x  = constant,  etc. 

tv*n»   T~D     =  shear  stresses 
xy   ap 

An,  En,  an,  bn =  complex  constants 

M  =  applied  bending  moment 

y0  =  half-height  of  opening  at  x  =0 

=  s  -  t  /  R 

I  =  moment  of  inertia  of  beam 

__=  bar  above,  indicates  conjugate  complex 

SCF =  stress  concentration  factor 

-  37  - 


I •■evelopiient  of  Th^or^tlc.-.!  Kolu'lon^ 

Consider  a  function:     z«  f  (  £  )  (1) 

which  transforms  continuously  points  in  the   £ (  a ,  p  }  plane 
into  points  in  the  z(x,y)  plane.  'Apt?  zm    x  /  iy  and 
K  =      e  a+iP  =   ew. 

In  particular,  se?  Fige.  XV(a)  and  (b),  points  on  the 
circle  of  radius  ea  ,  with  a  held  oonstant,  defined  by  a 
position  angle  p  ,  will  transform  into  points  (x,y)  on  a 
closed  curve  defined  by  the  functional  form  (1).   The  radial 
line  defined  by  p  =  constant  transforms  into  a  curve  orthogonal 
to  those  for  which  a  s  constant.   The  case  of  a  -   0  is  im- 
portant; it  defines  the  unit  circle,  f      ,    in  the   K   -plane 
and  transforms,  by  the  proper  selection  of  the  form  of  (1), 
into  the  shape  of  the  opening  to  be  studied.   A  special  symbol 

X   will  be  assigned,  to  those    ^-values  which  lie  on  this 
circle: 

X  =  elP  (2) 

so  that  the  conjugate  of  X  :  X"  ■  e*"^  P  ■  1/  X. 

The  transformation  defined  by  Equation  (1)  describes  the 
relation  between  an  increment  dz  in  the  z-plane  and  an  increment 
dw  in  the   ^-plane.   The  ratio,  expressed  in  polar  form: 

25  =  he  **  (3) 

dw 


3  -  Thin  development  follows  that  of  (B)  and  was  largely  guided 
by  1  r.  J.  S.  Brock. 
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defines  h  the  stretch  ratio,  unci  X  ,  the  an;  le  between  the 
tangent  to  the  curve     s constant  and  the  x-axis;   ape  [lo]  , 
Article  2.32.   h  and  %     are  obtained  by  differentiating  (1): 

dw   cL    dv   cL  <-, 
Equating  (3)  and  {■':): 

where  the  prime  indicates  differentiation  with  re?:ppct  to  the 
variable  in  the  parenthesis. 
Thus,  the  conjugate  of   (5)  Is: 

provided  (1)  has  real  coefficientp  and  parameters. 
Multiplication  of  (5)  and  (6)  gives: 

h2  „  t«(-4)«M  5)  5  I   =zM  S  )z»(  5  )e2  a      (7) 
The  division  of  (5)  "by  (6)  gives: 

In  general,  in  plane  problems  of  elasticity,  the  tech- 
nique which  follows  can  be  applied  to  any  real  rational  map- 
ping function  ?.    =f{  £  };  that  is,  where  21  is  the  ratio  of  two 
polynomials  in   £   havin  real  coefficients. 
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Consider  the  mapping  function  cur-  to  Greenspan  t7J  • 

»  -  8  *  /  |  /     ^  (9) 

where  s,  t,  and  a  are  real  constants  that  are  determined 

from  the  -proportions  of  the  opening,  and  ?,   =  x  /  iy  and 
?  b  e  a+ip  # 

Then:   x  /  iy  «  sea  •  e1?  /  te~a  •  e""1?  /  Re~-a  •  e~?1P 

m   (sea  /  te"a)cor  B/  Re~3acor3B  /  i(sea-  te"*a)sin  9 


-  iRe~3a  sin3  a- 

uatlng  the   real  and  iniaglnary  parts   separately: 

x  s   Uea  /  te~a)   coep     /  Re~^a  cos3   B 

(10) 

y  =  (sea  -  te~a)  sinB  -  Re~--a  ein3  b 
'hen  a  B  0,  Equations  (10)  define  the  general  ovalold  opening: 

X  =  (s  /  t)  COSo   /  R  0O83  R 

(10a) 
y  =  (s  -  t)  elnB  -  R  sin3  p 

Several  special  oaseo  are  of  interest: 

(a)  T'r.en  H  =  0,  an  ellipse  is  defined. 

(b)  when  R  b  t  =  0,  a  circle-  of  radius  n  is  defined. 

(c)  When  t  =  0  and  F  <  0,  a  scuare  with  fillets  at  the 
corners  is  approximated. 

(d)  When  t  =  0  and  R  ?•  C,  the  square  is  rotated  ^5°. 

In  the  general  case,  by  the  proper  choice  of  s,  t,  and  R,  the 
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equations  ca  be  made  to  approximate  the  common  structural 
openings ,  specifically  the  one  eora.oosed  of  a  rectangle  vith 
f eraiclrclec  at  each  en  . 
From  (9):    z»  (  £)  "  s  -  -«  -  2£_  (.1) 

Nov  h  and  %  can  be  determined  from  (7)  and  (8).  '.'hen 
a  =  0,  corresponding  to  the  opening,  the  subscript  zero 
is  used  ajvJ  from  (7): 

y&  -  «•<   X)«M   X)   «  (8  -  ^  -  2g)(B  -  *         2*   ) 
0  x^       x^  x^        xr 

(12) 

=  s2    +  t2    -*9R2  -    (2«t  -   6Rt)cos2    P  -   6Rflcos4P 

In  order  to  obtain  stresses   in  plane   problems  of 
elasticity  by  using   the  complex  variable  method,    two  potential 
functions  must  be  determined.      These   functions  can   b<*   expressed 

i    -    • 

$U)    -       £     ARzn     /         £    ,^ 

(13) 

vhere,  An,  an,  Fn,  and  bn  are  constants,  possibly  complex. 
7he  solution  Is  thus  re.dubed  to  the  evaluation  of  these  con- 
stants.   ih?>  stress  condition  at  on<=-  boun-: 'ary ,    infinity,  vill 
determine  A  and  B  ;  the  strorf  condition  at  the  other  boundary, 
the  opening,  vill  determine  a  and  b  .   Once  the  t-o  potential 
functions  ^re.  completely  determined,  the  stresses  in  the  plane 
may  be  determined  by  differentiation.    ub,  these  functions 
-.r-3  counterparts  of  the  Airy  strers  function.   In  fact,  the 
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Mry  stress  function  corresponding  to  these  potential  functions 
is  given  by: 

F(x,y)  =  Real  Part  of  T  z  $(z)  /  J"$(z)dz]    (14) 
The  stresses  can  be  found  in  terms  of  the  two  functions: 


ov  /  a  =  2f  §«(z)  /  |«(z)l 

x    y    *-  -* 

-     0X/    2iTxy    -<2[l     $«(Z)    /     Ji(E)J 


(15) 


The  usual    transformation  equations      . ive: 

(o{3    *   °a  +  2lTap)Q""i=  ay  *    °x  +   *lTxy  (l6> 


By  combining    Equations    (8),    (15),    and    (16): 

o      +  oQ      =2   £$•(*)   /    ImzJ)  (1?) 

where  the  functions   of   z  arf;   to  be  expressed   in   terras   of 
S  (a  ,  p  ): 

<§U)  =     $(*<*>]  =  <p<«>  (l8) 

$(z)    =      $§(£>]  -  <p  <«) 
Differentiation  ie   with   respect   tj   the  variable   In   parenthe- 
B  R  r  : 

U  -   net.   [ld\    ,   p.   135 


z'k) 

ft'(z)  =  SiiSi  (19) 

-'(5) 

If,  at  the  opening,  the  normal  stress  a   and  the  shear 

a 

stress       t   0   are  zero,   vhich  is   the   case  at  a   free  boundary, 
ap 

it   can  be   shorn  that: 

<p    (    X)   +   -SiiJL    ^    »(J.)      +   f    (-1-)    =   0  (20) 

The  functions  (p   (  £)  and  ^  {  £ )  may  be  evaluated 

rith  th~-  use  of  the  Cauohy  Theorems  regarding  poles  and  resi- 

5 

dues. 

From  [12]  ,  these  are: 

Theorem  I.   If  f{  £)  is  continuous  in  the  closed  region 
[  ^\   -  1  and  analytic  in  the  interior  vlth  the  possible  exce  >- 
tion  of  the  point  Z,   *  0,  where  f(  4)  has  the  form 

f(  K)   =  &  /  §§  /  •••■  /  |fi  /  gU  ) 

and   g(  £  )    is   analytic;    then: 

for     |£  |  >     1 
\<]\ere      y    is   the  unit  circle    y  =      e3-  ft 


-    K"    - 


Theorem  II.  Ii    f(£  )   is   continuous   in  the  cloned    region 

\£\^    1  and  analytic  in   th<=   r^c-ion  exterior     o      *•      ,   with 

the  possible  exception  of  the  point      ^  =  oo    ,    where   f (  K  ) 

has    the   form 


2    «  1  -/r  *-<,  0     12  n 

for  |^|  >    1 

Consider  Equations    (13): 

S...I    \  -       4       i      ,  .a     n-2     v        n(n+l)an 

Now  lpt  z  approach  Infinity   in  the  z-directlon: 

The  boundary  conditions   at  x      =    °°are: 

-  ^  - 


*.=  <*> 


K  =  °° 


(21) 


M 

x  I   '  u  y  '       l  xy 


Inserting  E  ; nations    (21)    Into  the   firet   of  Equations    (15) 
along  --ith  these  boundary  conditions   gives: 

-  £y  =    2[Z    ni^x  ^   iy)11""1  /     £nAn<x  -   iy)""^^ 

This  requires  that  A  n  =  0  for  all  values  of  n  except  2.   For 
n  =  2 


-  P    =  2  ^2iA2-   -  2iA2y] 


A~  =  1  JL 
2     81 

Inserting  Equations  (21)  into  the  second  of  Equations  (15) 
gives: 

|y=  2  {(-  -   ly).2.1-A2/  Z   nB„(x  /  iy)""1]  ^  _  „ 

This   requires    that  B      =  0   far  all  values  of  n   except  2.     For 
n    =  2 

|y=2[£y/    2lr>2v] 

Bf        _i«fy=_A2 

Therefore,  for  pure  bending,  regardless  of  the  shape  of  the 
opening,  the  potential  function?  are: 


$(Z)  =  *-y2/     Z  anA„ 


(22) 


$U>    --gar 2/    f    V2n 
On  the  boundary   of   this    opening: 
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81  A      7  '    (£  A+i+*^)n 

X      X  ^ 

X  L     '^  ^  '      (5  +t   X    +R  X     )n+1J  X 

^      (a)   =-  i£(B   X+  i   +£,)2+      f 


81  x     x^  '      r«x  *  5  JLi 


bn 


(f   +tX+HxV 


The  boundary  condition   (20)   bpcomes: 

4.    .,2 


<p    (a) , —  (-    +t    A+  R   A'  )♦! D ^~n+T 

A"3  Ul   X  '     (b    +u   +  R   x3)n+IJ 


n ,-» 


+        i*(§    +tX   +  RA    ■>)%   £ 

81  '     <j5   +tX   +  R  x3)n 

(23) 
Expanding  and   collecting  terms: 

*    (X)  -  -- (s2  +t2  +F2  -st)   -  i-(2st  +2Rt  -   t2  -   2Kb)   x' 
41  81 

i''h.p  orj-M    ^    I.1HT32   1^  lM.2st    +2Rt    -    s2 

-   --12^   -   2Rt)   X      +  --R     X      -   --( * ) 

?5l  81  ol  a^ 

IK         >Ra         y  nSD  .b   X4  +t  X2+   Rx 

>     (-  +t  x   +R  x-M 

At     £   =0,      cp  (   £  )   has    the    fon: 

,     i       >        D..-       •)      lM/2flt  +2Rt  -s2\    .  -Hs,  % 

v         ^  4      +oI  ^2  81  £**      F,^ 
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From  Theorem   I: 


T71    J     ~X™  81          ™2~                  31        Fr'        s2if 

r  4 

for   |s|>  1 

-Jw.        fxsixl2h  =  .   i2S(IfJi--2Rt-l£f )   -  £2        Sgs 

2*  1     °^        X-£  81               ^  ^                   81    *     ^ 

y  for  Kl  >    1 

At     £  =    oo     ,      <p  (£  )   beiB   the  form: 

«  T  :             «  0                   »      D       .^r 

81  '      £K 


I'rora  Th^or^m  II: 


i      r  slxi^x  - 

2*  i     J       xZZ 

-    <p<£   )   +  iS   e2      £2 
81 

Y 

-J—     C  xsihlzx  = 

2   *i     ^        X-4 
Y 

-   Sep     U  >    +  ~  a*"  S 

for     UI>1 


n 


J  for     UI>1 


From  the   firrt   of   the   expressions  derived  from  each  theorem: 

81  *-  ^2  ^  4J  e2^ 

From  the  second   of   the  expressions   derived    from  each    theorem: 
t      \        1Mr*2   x3      2et+  2Rt  -s2      2Rs"| 

C omparing    ( 2  k )  and   (25): 

a,    =  0 

-  if7  - 


Putting    T  =  t  /  S^  -  in 

r    2 


f  u ,  -  g«  [i£  ♦£  *A]  U6) 


But  the  conjugate  of  Equation  (20)  is: 
Ueini  Equation  (26),  this  becomes: 


?T        4FU  (%   +  tX    +T;   ^     ) 

41    1*2*  J       41  ^3      \5(E.i,  .  2g) 


-  125  ["---  +TA2+F  x'H  +    HH 
41     I  2  *  J        41 


+    y  (  X)    =  0 

,_  -.4  r.m  4R  _ 

,<x>    .  m  [  •      +TX2  +  R  X4  -  f*       if  ~  S-<s    +t   x2    ^     *H 

At  ^  =   0,    ^(    £)   has   the   form: 

From  Theorem  I: 


JL      /"  x 


'(--^  for    |«|>   1 

2*  i    ^         *<  41      2^ 

Y 

At   £  =    °°  ,  iM  £  )  has  the   form: 

0  v  4     0«       4B 
/     x        lKsr,     2      „     4        Sc»       ~ci    ~  "^    ,         ^y2        ^    >J 
V  (£)  = T^     +  R   £ 5 •-- A-    (s    +tC    +  R   £  ) 

41  u  S£     -t£     -   3R  -J 

From  Theorem  II: 


2    TTi     v/ 


s£4   -    2T   -    4R 
*(X)dX         .    ,     ^      1KB      m  2        _      4  ^  K2 

X-£  ^1  s£      -    t£  -    3R 

Y  (s  +  t£2  -t-  R£4)J         for    \^|>    1 
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Therefore: 


<p  («)  .  £&[-■  +  --  :+p^ 


*  *« 


■    -2      -   -J 


a 


£'   -t*     -    3! 


^-(8+tC+} 


i«*>| 


(27) 


lining  the  first   of  Equations    (17)    ^nv   no tin*    that  at   a   free 
boundary    a  =       0,    oq   =     0 

P    a=0  ^  a=0         Lz'(\)      z'(f-)J 


But 


,.,  s      IMS  TsA2      T        R~l  -      -  tl\ 

<*>(x)=i  zii"  IT**  ~-  +  "nJ   a       <p(^)=- 

<p«(X)  =  lHfifBX  -  2JC  -   iifl  and    «   »(!)  = 


z'(X)    =e   -  ^ 


X" 

V 


and      i»(|) 


■  5rt2x2+  iX  +"xj 

.  iMfiffi   -2TX3-W] 

2,1    L  X  J 

t  s ■-   t   X  -   3H  X' 


Thus : 
( 


P    a=0     L"1  B      .   .t         2g  41  g  „   tx2  .    3Rx4     J 

X  X' 

Pa=0        x    L(P    +t   +9R    )   +(6Rt  -   2st)cas2p  -   6RBCos4p  J 


(23) 


(0  )    Is  a  maximum  vhen  p  =  90°. 
P  a=0 


Therefore: 


SCF  = 


(oR) 


c 


-  &fl 


_   6 

"  y0 


2e  -3R8-12R^+2t(4R-t)(l+g) 


(s  +t  -  3R)' 


(28a) 
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I  f  ve  noi  ■  s  u ;,£-.  ti t ut e  Eq  ua  isl one  ( 1 ) ,  ( 9 ) ,  (26),  (:•:?),  and 
(11)  into  Equation  (14),  integrate  and  combine  terras,  vp 
get   the  Airy  streee   function: 


(afp)   =  Kfi  |f  (62«Hrt-2Rt+3K2)ea 

4i  Lu 


+(3rT  -  fi-  -2T*  -4R^)e 


_c    +3a_^        -3a  -5  a.      ~    -7a  "I 

-  £-e        *fte       -RTe        +  P/e        J    ninp 

t,st  +  Rt  -sTx    5a      ,st        2sT        4nt,    -3a.   ^    -a 
( . )e     +C7--5 5-)e       +*Te 

st    a  -bal  r        -5a     Re    -a     R      -6al  ^ 

-  — e    +  Rte       J  sin3  p   +pse       -  -~e     -  — e       J  ein5  p  V 
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Evaluation  of  Parameters 

It  is  necessary  to  of-ter-nine  the  values  of  R,  s,  and  t 
that  will  ensure  that  the  mapplnc  function  (;)  will  suit  the 
actual  opening  vhen  gC«  C.   Several  conditions  may  be  im- 
posed on  (10a)  in  order  to  determine  B,  s,  -.in-     t  Such  *s: 

(a)  s  /  t  /  R  «  &  horizontal  dimension  ■  >a 

(b)  s  -  t  •/  R .  «  ;-.  vertical  (linens ion  *  -;b 

(c)  Area  of  opening  ■  (a  -  b)b  /  ..bntt  **JydiX 

(d)  Radius  of  curvature  at  end  ■  "b 

^ouationR  associating  the  aesired  variables  with  the 
geometric  properties  of  the  opening  must  nov  be  developed  to 
provide  a  eyeteraatio  procedure  for  evaluation.   To  simplify 
this  orobl^n  we  introduce  the  Greenspan  notation: 

(30) 

*  (P  -  Ci)  a  t 

Or: 

8  /    t  *  P 

(30a) 

s  -  t  »  fl 

Then,    Equations    (10a)   become: 

x  •   F  oosr   r   R  cos    3b 

P  (10b) 

y  «  C£  cinp   -   R  sin  3? 

Imposing   the  condition  of  equal  ar°as,    T'e  obtain: 

h    J     y  dx  »   (a  -   b)b  /  £b2*      for  a  >    b  (31) 

o 

Substituting  (10b)  into  (3D,  integrating,  and  simplifying 
gives: 
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%     (p     .   3p2)    ,  (a  .    L)b  /  .^2^  (  ?2) 


Reference   to  any  text   on  differential   calculufl    shove   that 
th?3  radius   of  curvature  at  any  point   is  given  by   the  ex- 
pression: 

Radius** H3) 

■««  *■  -Si  =iv&  anfiy"=  fi1  -IV -Si 

But ,   -^  =  Ci  COS  R  -  3R  cos  3  8 
dp        r  r 

and   flJ- —  1 

"      P  sin  e /  3H  sin  ?  b 

Hence:      y«  -    -  SLS9UL^2LSSSL3J  (34) 

?  sinP  /  3H  sin  3  P 

„  -PQ  /  2?R*-3R<F  /  cjcos   4P/  6R(P-;)cos   2  0 

and  y    =    4  •— — -.--*.- -.- —_——•. --— . 

(P  sin  b  /  3R  sin  3  p  ) 

(35) 
Substituting    (3^)   and   (35)   into   (33)   and    simplifying  gives: 

Radius  -  iili^lliizhfz{2ldz2S^i^ll^2iZh23i?t<^l22^l!tM2^l. 

-?c,/27H2-3H(P/-Joos4B/6n(P-c;)cos2P 

Now  b=  0  for  the  longitudinal  end  of  the  opening. 

Therefore: 

Padlus  =i-~--21}2-  =  fb  (37) 

P  /  9H 

The  remaining  two  equations  are  b imply: 


p 


/  n  =  ^a  (33) 


Q  /  R  =  iVb  -(39) 


5^  - 


isidering  only  Equations  (32),  (37),  (38),  unr    (39)  as 
a  system  of  equations  for  the  determination  of  P,  %,  and  R, 
we  now  choose  b  =  ?  and  involve  the  ratio  o:  length  to  depth 
of  opening,  a/b    .   The  system  of  equations  become: 
I?  /   R  =  a/b 

;  /  n  =1 

P.  -  31' "  =  -(a/b  -  1)  /  1 

ia-=j5BJ:  =  1 

P  /  9H 

Any  system  of  iaore  equations  than  there  ..ire  unknowns  is  In- 
determinate.  Therefore,  the  four  possible  combinations  of 
(kn)  taken  three  at  a  time  were  solved  for  P,  Q,,    and  R  for 
the  particular  case  of  a/b  =  Z.      This  particular  case  was 
chosen  since  it  had.  been  used  by  Greenspan  t  7  J  as  an  il- 
lustration and  a  ready  comparison  vas  available.   The  outline 
of  the  opening  as  defined  by  there  four  solutions  are  plotted 
in  Fig.  XVI  and  compared  with  the  actual  opening;  for  vhioh 
the  fit  was  desired* 

Because  the  ratio  of  depth  of  opening  to  depth  of  be^n 
is  an  important  criterion  in  the  determination  of  the  strese 
concentration  factor,  the  depth  of  opening  seemed  the  post 
logical  first  choice  to  be  <=>ld  constant.   Next  in  importance 
as  governing  variable  Was  the  length  of  opening.   To  complete 
the  system  of  equations  for  solution  for  P,  Q,  and  R  either 
the  radius  of  curvature  at  the  longitudinal  end   or  the  area 
must  be  chosen.  A  comparison  of  the  fit  showed  that  either 
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would  be  satisfactory,  and  a  comparison  of  the  two  values  of 
stress  concentration  factor  examined  with  Greenspan's  value 
as  the  yardstick  shoved  that  either  v  >uld  be  satisfactory. 
A  good  case  could  dp  sriade  for  holding  the  pnd  radius  of 
curvature  equal  to  the  half  height  of  opening  since  this 
system  gave  valuer,  of  stress  concentration  factor  that  were 
slightly  high  -  pessimistic.   When  plotted  on  a  curve  similar 
to  Fig.  Ill,  ho- -ever,  the  closeness  of  fit  of  values  calcu- 
lated by  holding  the  area  constant  was  slightly  superior  to 
those  values  calculated  by  holding  the  end   radius  of  curvature 
equal  to  the  half  height  of  opening.   ?he  --.re  of  calculating 
the  valuer  of  P,  \.,  and  P  -hen  holding  the  area  constant  was 
a  further  cone id era tl on. 

Consequently  a  system  of  equations  consisting  of  the 
first  three  of  Equations  (&0)  vas  chosen: 

P  =  a/b  -  R 

=  1  -  ':  (4C) 

PQ-3R2  =  -(a/u  -  1)  /  1 

Frora  which: 

R  =  -  .(a/b/1)  /  i   /(a/b/1)2-  3<V*-  l)U/b-l)    (40a) 

These  valuee  of  P,  Ci,  and  P.  are  plotted  versus  a/b  in  Fig. 
XVII. 

This  should  not  be  construed,  to  mean  that  th^se  values  are 
precise;   the:'  are  merely  the  best  choice  that  the  authors  could 
make  with  existing  information.  Somp  other  geometric  property, 
such  as  the  peripheral  length,  may  warrant  consideration. 
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